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Abstract. Hard and soft QCD dynamics are both important in charmonium hadroproduction, as 
04 | presented here through a next-to-leading order QCD matrix element calculation combined with the 

■ colour evaporation model. Observed x F and p ± distributions of //y/ in hadroproduction are repro- 

duced. Quite similar results can also be obtained with a Monte Carlo event generator where cc pairs 
, are instead produced through leading order matrix elements and the parton shower approximation 

of higher order processes. The soft dynamics may alternatively be described by the soft colour in- 
teraction model. We also discuss the relative rates of different charmonium states and introduce an 
improved model for mapping the continuous cc mass spectrum on the physical charmonium reso- 
i— i ■ nances. 

> : 

I/-) ! The theoretical description of charmonium production separates the hard and soft 

parts of the process based on the factorisation theorem in QCD. Thus, we first consider 
the perturbative production of a cc pair at the parton level and then the non-perturbative 
(Sj , formation of a bound charmonium state [jl]]. 

Perturbative QCD (pQCD) should be applicable for cc production, since the charm 
^ ■ quark mass m c is large enough to make (X s (m%.) a small expansion parameter. The 

leading order (LO) processes are gg — > cc and qq — > cc. The next- to-leading order (NLO) 
Qh| processes, i.e. @(a%), include the emission of a third parton and virtual corrections 

(where divergences are properly cancelled). The full NLO matrix elements, with explicit 
charm quark mass, are available in a computer program [Q] giving total and differential 
cross sections. 

An alternative description of the pQCD production of cc pairs is given by the PYTHIA 
Monte Carlo, where all LO QCD 2^2 processes are included with their correspond- 
ing matrix elements and the incoming and outgoing partons may branch as described by 
the DGLAP equations. A cc pair can then be produced as described by the LO matrix 
elements for qq — > cc and gg — ► cc (with explicit m c dependence) or in a gluon splitting 
g — > cc in the parton shower. 

The main free parameter is the charm quark mass m c , taken as m c = 1.5 GeV in the 
NLO program and m c = 1.35 GeV in Pythia. In both approaches, the factorization and 
renormalization scales are taken as the average transverse mass of the c and c. 

The formation of bound hadron states occurs through processes with small momen- 
tum transfers such that a s is large and prevents the use of perturbation theory. The lack 
of an appropriate method to calculate non-perturbative processes, forces us to use phe- 
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nomenological models to describe the formation of charmonium states from perturba- 
tively produced cc pairs. The Color Evaporation Model (CEM) [Q] and the Soft Colour 
Interaction (SCI) model |6|] are based on a similar phenomenological approach, where 
soft colour interactions can change the colour state of a cc pair from an octet to a singlet. 
They employ the same hard pQCD processes to produce a cc pair regardless of its spin 
state. A colour singlet cc pair with an invariant mass below the threshold for open charm 
(m cd < 2m D ) will then form a charmonium state. 

In CEM [Q, 0, ||, ^|] the exchange of soft gluons is assumed to give a randomisation of 
the colour state. This implies a probability 1/9 that a cc pair is in a colour singlet state 
and produces charmonium if its mass is below the threshold for open charm production, 
m cd < 2m D . The fraction of a specific charmonium state /, relative to all charmonia, is 
given by a non-perturbative parameter p- (Pj/y = 0.4 — 0.5) [Q]. 

In SCI [§, ||, [K]] it is assumed that colour-anticolour, corresponding to non- 
perturbative gluons, can be exchanged between partons emerging from a hard scattering 
and hadron remnants. The unknown probability to exchange a soft gluon between parton 
pairs is given by a phenomenological parameter R. These colour exchanges lead to dif- 
ferent topologies of the confining colour string-fields and thereby to different hadronic 
final states after hadronisation. The mapping of cc pairs, with mass below the threshold 

for open charm production, is here made based on spin statistics resulting in a fraction 

r 

of a specific quarkonium state i with total angular momentum /■ given by f i = 

where T = (27-+ including a suppression of radially excited states through the 
main quantum number n-. This model was found to give a correct description of the 
different heavy quarkonium states observed at the Tevatron []5p. 

The complete models are formed by adding the CEM or SCI models for the soft 
processes to any of the descriptions for the hard pQCD processes. The first model we 
label CEM-NLO and is the combination of the CEM model with the NLO program. The 
second model is CEM-PYTHIA, where CEM has been implemented in Pythia version 
5.7 [J3p. The third model, SCI-PYTHIA, is to use the SCI model as implemented in 
Pythia 5.7. Further ingredients are the intrinsic k,, due to the Fermi motion of partons 
inside the initial state hadrons, and soft p T in soft gluon exchange that neutralize color. 
Both effects are modelled by a gaussian distribution of width 0.6 — 0.8 GeV used in 
Pythia and in the NLO program. 

Comparing these three models we can separate different effects. With CEM imple- 
mented in the NLO program and in Pythia, we can compare the pQCD contributions, 
namely NLO versus LO plus the parton shower approximation of higher orders. Hav- 
ing SCI and CEM implemented in Pythia, we can explicitly compare these two non- 
perturbative models and see to what extent they can account for observed soft effects. 

Detailed comparisons between the models have been done as well as extensive com- 
parison with data, both from fixed target experiments and the Tevatron collider [jlp. Here 
we limit ourselves to proton beams. The targets are different nuclei, but the experimental 
results are rescaled to the cross section per nucleon. Thus we compare directly with our 
models which do not include any nuclear effects but treat hadron-nucleon interactions. 

Fig. H] shows x F and p± distributions of the produced J/y for proton beams of 
different energies. As can be seen, the data are approximately reproduced, both in shape 
and normalization, by all three models. Looking into the details of the x F distributions, 



one can observe that the model curves fall less steeply than the data and therefore 
overshoot somewhat at large x F . The observed p ± distribution is better reproduced, with 
only small differences between the models. 




FIGURE 1. Distributions in ir and p\ of J /iff produced with proton beams of energies 800, 530 and 
300 GeV on fixed target. Data Q [11 |3j compared to CEM based on NLO pQCD matrix elements, 
and CEM and SCI based on LO matrix elements plus parton showers in the Pythia Monte Carlo. 

Having CEM combined with different treatments of the pQCD production of cc, we 
can now investigate pQCD effects in more detail. Fig. ^| illustrates this for the case of 
800 GeV proton energy, similar conclusions can also be drawn for other energies and 
beam particles. For the x F distribution in Fig. ^a, the full NLO result and that based 
on LO+PS agree reasonably well. The NLO corrections are very important, as we see 
by comparing the LO and the full NLO results. In the LO+PS result, however, the PS 
contribution is unimportant for the overall cross section which is dominated by the LO 
cc production. The agreement with the NLO result is here obtained by using a lower 
charm mass, m c = 1.35 GeV. We have cross-checked this within the NLO program, 
where the full result is essentially reproduced by the LO part if this lower mass value is 
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FIGURE 2. Distributions in x F and p\ of J/\jf (in 800 GeV proton on proton as in Fig. [I]) from 
variations of the pQCD treatment. CEM based on the NLO program with m c = 1.5 GeV: NLO and LO 
matrix elements, NLO with no intrinsic £ ± . CEM based on Pythia with m c = 1.35 GeV: LO matrix 
elements plus parton showers (PS) and PS contribution shown separately. 



used. This demonstrates that the NLO correction is essentially an overall ^-factor from 
soft and virtual corrections. For the p ± distributions in Fig. |]b, the NLO program gives 
a p, distribution with a much larger tail at large p ^ , but it is still substantially affected 
by the inclusion of the intrinsic k ± at the limited values of p j_ accessible at fixed target 
energies. The p, distribution resulting from the LO+PS in the PYTHIA approach, is at 
high-p ± dominated by cc from gluon splittings in the partons showers, whereas the bulk 
of the cross section comes from the low-p ± region where the LO diagrams dominate. 
The total LO+PS result, which also includes a gaussian intrinsic k ± with the same width 
0.6 GeV, agrees quite well with the NLO result. 

Data on \j/ production provide an additional testing ground for the models, which 
produce all charmonium states with the same dynamics. A comparison made in [|l|] 
shows that all models account quite well for the shape of the distributions. The proper 
normalization of CEM is obtained by chosing p^, = 0.066. The spin statistics used in 

SCI predicts only a factor two suppression of y/, and must be lowered by an additional 
factor four in order to reproduce the data. This has prompted us to develop a more 
elaborate model for turning cc pairs into different charmonium resonances [[Ip, which is 
briefly described here. 

The cc pair is produced in a pQCD process with a continuous distribution of its 
invariant mass m cd and must be mapped onto the discrete spectrum of charmonium 
states. The soft interactions that turn the pair into a colour singlet and form the state, 
may very well change its mass by a few hundred MeV, which is the typical scale of the 
soft interactions. We model this by a gaussian smearing of a few hundred MeV. The 
probability to end up in a specific resonance, shown in Fig. [3|, is then proportional to the 
superposition of this gaussian with the resonance peak, times the corresponding spin- 
statistics factor. The smearing of m cd across the threshold 2m D for open charm, implies 
non-zero contributions for charmonium also above the DD threshold as well as some 
open charm production for m cd originally below this threshold. 




FIGURE 3. Probability distributions for the different charmonium states as obtained in the model with 
gaussian smearing (o sme — 400 MeV). The resulting total probability for charmonium production and the 
remainder as open charm production are also shown. 



By folding these charmonium probability functions with the distribution in m cg ob- 
tained from pQCD, one gets the cross section for a given charmonium state. Applying 
this mapping procedure to the CEM model we obtain the results in Fig. |j. As opposed 
to the simple spin statistics factor, this model gives a reasonable description of the ob- 
served ratio of i// to J/xy production and fractions of J/yr produced directly, coming 
from decays of % c states and from i//. In particular, the model gives a characteristic 
energy dependence of the kind indicated by the data. 

In summary, both hard and soft QCD dynamics play important roles in the production 
of charmonium states in hadronic interactions. The cc pair production in pQCD have 
substantial higher order contributions, with a factor two increase of the total cross section 
from NLO corrections. These come mainly from soft and collinear gluon emissions 
combined with virtual corrections and can be effectively accounted for by an overall 
^-factor. This supports to the use of the Pythia Monte Carlo with LO matrix elements 
and a reduced charm quark mass to increase the cross section correspondingly. The high 
/? ± tail of the cross section is, however, dominated by higher order tree diagrams in the 
NLO matrix elements and in the parton showers of the Monte Carlo approach. 

The non-perturbative formation of J /y, can be described by the Colour Evaporation 
Model and the Soft Colour Interaction model, where cc pairs in a colour octet state 
can turn into a colour singlet state by soft gluon exchange. A simple spin statistics 
factor is not sufficient for a proper description of other charmonium states, but our more 
elaborated model to map cc pairs onto the physical charmonium states improves this 
situation. 

To conclude, the main features of hadroproduction of charmonium can be described 
in these models combining pQCD and effects of soft colour exchanges. This shows, in 
particular, that these models for the soft QCD dynamics contain the essential effects and 
therefore improve our understanding of non-perturbative QCD. 
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FIGURE 4. The ratio of y/ to J / y/ production (times their branching ratios for decay into ju + (left) 
and fractions of J / y/ produced directly, and coming from the decay of % c and y/' states (right) in hadron- 
hadron interactions of cms energy y/s. Data [|l5| |l6| TTl |l8| compared to simple spin statistics and to our 
model with different gaussian smearing widths applied to CEM. 
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